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We study the effect of a magnetic insulator (Yttrium Iron Garnet - YIG ) substrate on the spin 
transport properties of NigoFe 2 o/Al nonlocal spin valve (NLSV) devices. The NLSV signal on the 
YIG substrate is about 2 to 3 times lower than that on a non magnetic Si02 substrate, indicating 
that a significant fraction of the spin-current is absorbed at the Al/YIG interface. By measuring the 
NLSV signal for varying injector-to-detector distance and using a three dimensional spin-transport 
model that takes spin current absorption at the Al/YIG interface into account we obtain an effective 
spin-mixing conductance G^_i ~ 5 —8 x 10^^ We also observe a small but clear modulation 

of the NLSV signal when rotating the YIG magnetization direction with respect to the fixed spin 
polarization of the spin accumulation in the Al. Spin relaxation due to thermal magnons or roughness 
of the YIG surface may be responsible for the observed small modulation of the NLSV signal. 


The coupled transport of spin, charge and heat in non¬ 
magnetic (N) metals deposited on the magnetic insulator 
Y3Fe50i2 (YIG) has led to new spin caloritronic device 
concepts such as thermally driven spin currents, the gen¬ 
eration of spin angular momentum via the spin Seebeck 
effect (SSE) [1], spin pumping from YIG to metals [2], 
spin-orbit coupling (SOG) induced magnetoresistance ef¬ 
fects [3, 4] and the spin Peltier effect, i.e., the inverse 
of the SSE that describes cooling/heating by spin cur¬ 
rents [5] . In these spin caloritronic phenomena, the spin¬ 
mixing conductance of the N/YIG interface controls 
the transfer of spins from the conduction electrons in N 
to the magnetic excitations (magnons) in the YIG, or 
vice versa [6-10]. The interconversion of spin current to 
a voltage employs the (inverse) spin Hall effect in heavy- 
metals such as Pt or Pd. The possible presence of prox¬ 
imity induced magnetism in these metals is reported to 
introduce spurious magnetothermoelectric effects [11, 12] 
or enhance G^^ [7]. Owing to the short spin-diffusion 
length A in these large SOC metals, the applicability of 
the diffusive spin-transport model is also questionable. 
Experimental measurements that alleviate these concerns 
are however scarce and hence are highly required. 

In this article, we investigate the interaction of spin 
current (in the absence of a charge current) with the 
YIG magnetization using the NLSV geometry [13-15]. 
Using a metal with low SOC and long spin-diffusion 
length allows to treat our experiment using the diffu¬ 
sive spin-transport model. We find that the NLSV signal 
on the YIG substrate is two to three times lower than 
that on the Si02 substrate, indicating significant spin- 
current absorption at the Al/YIG interface. By vary¬ 
ing the angle between the induced spin accumulation 
and the YIG magnetization direction we observe a small 
but clear modulation of the NLSV signal. We also find 
that modifying the quality of the Al/YIG interface, us¬ 


ing different thin-film deposition methods [4], influences 
G^^ and hence the size of the spin current flowing at 
the Al/YIG interface. Recently, a low-temperature mea¬ 
surements of a similar effect was reported by Villamor et 
aL[16] in Co/Cu devices where G^^ ^ was 

estimated, two orders of magnitude lower than in the lit¬ 
erature [4, 8]. Here, we present a room-temperature spin- 
transport study in transparent Ni8oFe2o (Py)/Al NLSV 
devices. 

Figure 1 depicts the concept of our experiment. A non¬ 
magnetic metal (green) deposited on the YIG connects 
the two in-plane polarized ferromagnetic metals Fi and 
F 2 , which are used for injecting and detecting spin cur¬ 
rents, respectively. A charge current through the Fi/Al 
interface induces a spin accumulation /i.s(f^ = (0, /is, 0)^ 
that is polarized along the y direction, parallel to the 
magnetization direction of Fi. This non-equilibrium /ts, 
the difference between the electrochemical potentials for 
spin up and spin down electrons, diffuses to both +x 
and —X directions of Fi/Al interface with an exponen¬ 
tial decay characterized by the spin diffusion length Aat. 
Spins arriving at the detecting F2/AI interface give rise 
to a nonlocal voltage Vn\ that is a function of the rela¬ 
tive magnetic configuration of Fi and F 2 , being minimum 
(maximum) when Fi and F 2 are parallel (antiparallel) to 
each other. 

For NLSV devices on a Si02 substrate, spin relaxation 
proceeds via electron scattering with phonons, impuri¬ 
ties or defects present in the spin transport channel, also 
known as the Elliot-Yafet (EY) mechanism. The situa¬ 
tion is different for a NLSV on the magnetic YIG sub¬ 
strate where additional spin relaxation due to thermal 
magnons in the YIG and/or interfacial spin orbit cou¬ 
pling can be mediated by direct spin-flip scattering or 
spin-precession. Depending on the magnetization direc¬ 
tion m of the YIG with respect to fig spins incident at 
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FIG. 1. (Color online) Concept of the experiment for m || /is. 
(a) A charge current through the Fi/Al interface creates a 
spin accumulation fis in the Al. The diffusion of fia to 
the F 2 /AI interface is affected by spin-flip relaxation at the 
Al/YIG interface. Scattering of a spin up electron {s=h/2) 
into spin down electron {s=—h/2) is accompanied by magnon 
emission {s=h) creating a spin current that is minimum (max¬ 
imum) when fia is parallel (perpendicular) to the magnetiza¬ 
tion of the YIG. (b) Profile of fia along the Al strip on a Si 02 
(red) and YIG (blue) substrate. The spin accumulation at 
the F 2 /AI is lower for the YIG substrate compared to that on 
Si02. 

the Al/YIG surface are absorbed (m _L /i^) or reflected 
(to II fig) thereby causing a spin current density Js(r/ 
through the Al/YIG interface [9] 

j/(to)|^=o = Grmxim x ^is) + Gi{mx ^ls)+Gs^ls. (1) 

Here to = {nix, niy, 0)'^ is a unit vector parallel to the 
in-plane magnetization of the YIG, Gr (Gi) is the real 
(imaginary) part of the spin-mixing conductance per unit 
area and Gs is a spin-sink conductance that can be in¬ 
terpreted as an effective spin-mixing conductance that 
quantifies spin-absorption (flip) effects that is indepen¬ 
dent of the angle between to and fig- 

When m \\ fig some of the spins incident on the YIG 
are reflected back into the Al while some fraction is ab¬ 
sorbed by the YIG. The absorption of the spin-current 
in this collinear case is governed by a spin-sinking effect 
either due to (i) the thermal excitation of the YIG mag¬ 
netization (thermal magnons) or (ii) spin-flip processes 
due to interface spin orbit effects or magnetic impurities 
present at the interface. This process can be character¬ 
ized by an effective spin-mixing interface conductance Gg 
which, at room temperature, is about 20% of Gr [5]. Be¬ 
cause of this additional spin-flip scattering, the maximum 
NLSV signal on the YIG substrate should also be smaller 
than that on the Si 02 . When rh J- fig spins arriving at 
the Al/YIG interface are absorbed. In this case all three 


terms in Eq. (1) contribute to a maximum flow of spin 
current through the interface. The nonlocal voltage mea¬ 
sured at F 2 is hence a function of the angle between m 
and fig and should reflect the symmetry of Eq. 1. 

Fig. 2(a) shows the scanning electron microscope im¬ 
age of the studied NLSV device that was prepared on a 
200-nm thick single-crystal YIG, having very low coercive 
field [2, 4, 17], grown by liquid phase epitaxy on a 500 
/rm thick (111) Gd 3 Ga 50 i 2 (GGG) substrate. It consists 
of two 20-nm thick Ni 8 oFe 2 o (Py) wires connected by a 
130-nm thick Al cross. A 5 nm-thick Ti buffer layer was 
inserted underneath the Py to suppress direct exchange 
coupling between the Py and YIG. We studied two types 
of devices, hereafter named Type-A and Type-B devices. 
In Type-A devices (4 devices), prior to the deposition of 
the Al (by electron beam evaporation), Ar ion milling 
of the Py surface was performed to ensure a transparent 
Py/Al interface. This process, however, introduces un¬ 
avoidable milling of the YIG surface thereby introducing 
disordered Al/YIG interface with lower G^i [18]. To cir¬ 
cumvent this problem, in Type-B devices (2 devices), we 
first deposit a 20 nm-thick Al strip (by DG sputtering) 
between the injector and detector Py wires. Sputtering 
is reported to yield a better interface [4]. Next, after Ar 
ion milling of the Py and sputtered-Al surfaces, a 130 
nm-thick Al layer was deposited using e-beam evapora¬ 
tion. Similar devices prepared on Si02 substrate were 
also investigated. All measurements were performed at 
room temperature using standard low frequency lock-in 
measurements. 

The NLSV resistance i?ni = Vni/I as a function of 
the applied in-plane magnetic field (along y) is shown in 
Fig. 2(b), both for Si02 (red and orange) and YIG (blue) 
samples. Note that the magnetizations of the injector, 
detector and YIG are all collinear and hence no initial 
transverse spin component is present. The spin valve sig¬ 
nal, defined as the difference between the parallel Rp and 
anti-parallel Rap resistance values, Rsv = Rp — Rap 
on the YIG substrate is about two to three times smaller 
than that on the Si02 substrate. This reduction in the 
NLSV signal indicates the presence of an additional spin- 
relaxation process even for to || fig. Assuming an iden¬ 
tical spin injection efficiency in both devices, this means 
that spin relaxation in the Al on the YIG substrate occurs 
on an effectively shorter spin relaxation length A at. To 
properly extract Xpf we performed several measurements 
for varying distance between the Py wires, as shown in 
Figure 2(c) both on Si02 (red diamond) and YIG (blue 
square) substrates. Also shown are dashed-line fits us¬ 
ing the expression for the nonlocal spin valve signal Rsv 
obtained from a one-dimensional spin transport theory 
given by [14] 

^ ^ _ , 2 ) 

(lv + l)[|^sinh(d/2AAr) -b cosh(d/2AAr)]' 
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FIG. 2. (Color online) (a) Scanning electron microscopy image of the measured Type-A device. Two Py wires (indicated 
by green arrows) are connected by an A1 cross. A charge current I from contact 1 to 2 creates a spin accumulation at the 
Fi/Al interface that is detected as a nonlocal spin voltage Pni using contacts 3 and 4. (b) The NLSV resistance Rni = Vni/I 
for representative YIG (blue) and Si02 (red and orange) NLSV samples. For comparison, a constant background resistance 
has been subtracted from each measurement, (c) Dependence of the NLSV signal on the spacing d between the injecting 
and detecting ferromagnetic wires together with calculated spin signal values using a ID (dashed lines) and 3D (solid lines) 
spin-transport model. For each distance d between the injector and detector several devices were measured, with the error bars 
indicating the spread in the measured signal. 


Here Rp = (1 — and ^ are spin area re¬ 

sistance of the ferromagnetic (F) and non-magnetic (N) 
metals, respectively. Xn and Xp are the corresponding 
spin diffusion lengths, ap {cr^) is the electrical conduc¬ 
tivity of the F (N), Of is the spin polarization of F and d 
is the distance between the injecting and detecting ferro¬ 
magnetic electrodes. Fitting the Si02 data using Eq. (2), 
we extract ap =0.32 and AAr^iOa =320 nm, which are 
both in good agreement with reported values [13-15]. 
A similar fitting procedure for the YIG data, assuming 
an identical spin injection efficiency, yields an effectively 
shorter spin-diffusion length Ajv yig= 190 am due to the 
additional spin-flip scattering at the Al/YIG interface. 
This value of Aat^yig therefore contains important infor¬ 
mation regarding an effective spin-mixing conductance 
Gs that can be attributed to the interaction of spins with 
thermal magnons in the YIG. When spin precession, due 
to the applied external field as well as the effective field 
due to Gi is disregarded, we can now estimate Gs by 
relating Xn,yig to AAr^siOa via Gs as (see Supplemental 
Material [19], Sec. I): 


\2 

^Af.YIG 


'^YSiOa 


^ A?’ 


(3) 


with A“^ = 2Gs/tAicrN [19]- Using the extracted values 
from the fit, (TAr=2x 10^ S/m and t2ii=130 nm, we extract 
Gs — 2.5 X 10^^ 12“^ m“^, which is about 25% of the 
maximum Gr ~ reported for Pt/YIG [4, 7] 

and Au/YIG [8] interfaces. 

To quantify our results we performed three- 
dimensional finite element simulations using COMSOL 
Multiphysics (3D-FEM) [19, 20] that uses a set of 


equations that are equivalent to the continuous random 
matrix theory in 3 dimensions (CRMT3D) [21]. The 
charge current j“(r) and spin current j'^{r), (where 
a S x,y,z), are linked to their corresponding driving 
forces via the electrical conductivity as 


^ apa\fVy.c\ 


(4) 


where yLc = + M 4)/2 and /Xs = (mt ~ are the 

charge and spin accumulation chemical potentials, re¬ 
spectively. We supplement Eq. (4) by the conservation 
laws for charge (V • j“(r) = 0) and spin current (V ■ jg = 
(1 —a|.)cr [/Xs/A^ -b wl X Ms]) where ujp = gfipB/fi with 
g = 2 is the Larmor precession frequency due to spin 
precession in an in-plane magnetic field B = (i?x, By, 0)^ 
and gs is the Bohr magneton (see Supplemental Mate¬ 
rial [19], Sec. II). To include spin-mixing at the Al/YIG 
interface we impose continuity of the spin current js at 
the interface using Eq. (1). The input material param¬ 
eters such as cr, A and ap are taken from Refs. 22 and 
23. 

The calculated spin signals obtained from our 3D-FEM 
are shown in Fig. 2(c) for samples on Si02 (red solid line) 
and YIG (blue solid line) substrates. By matching the 
experimentally measured NLSV signal on the Si02 sub¬ 
strate with the calculated values in the model we obtain 
ap = 0.3 and Xn =350 nm. Using these two values and 
setting Gs — 5 X 10^^ well reproduces the mea¬ 

sured spin signal on the YIG substrate. This value of Gs 
obtained here is consistent with that extracted from our 
ID analysis based on Eq. 2. Hence, the interaction of 
spins with the YIG magnetization, as modeled here, can 
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FIG. 3. (Color online) (a) Nonlocal spin valve resistance Rni of a Type-B device with d=500 nm between injecting and detecting 
Py wires and tAi=130 nm. A constant background resistance of 117 mfl was subtracted from the original data, (b) Angular 
dependence of the NLSV signal in the parallel and antiparallel configurations. The AP curve is average of 10 measurements 
and that of the P state is a single scan. Both resistance states exhibit a cos(2a) dependence on the angle between m and fis- 
The black solid lines are calculated using the 3D-FEM model for Gr = 1 x that show a percentage modulation of 

only 12% corresponding to the green curve in (c) SRsv /Rsv is plotted. The angular dependent measurement in (b) is from 
a device for which complete set of measruements were peformed. A spin valves measurement as in (a) was also performed for 
another device with d = 300 nm. 


capture the concept of spin-mixing conductance being 
responsible for the observed reduction in the spin signal. 

In the following we investigate the dependence of Rn\ 
on the angle a between fig and m. We rotate the sam¬ 
ple under the application of a very low in-plane mag¬ 
netic field B <5 mT, enough to saturate the low-coercive 
(< 0.5 mT) YIG magnetization [4, 5] but smaller than 
the coercive fields of Fi and F 2 (~20 mT). This con¬ 
dition is important to maintain fixed polarization axes 
of /is, along the magnetization direction of the injecting 
ferromagnet, and also have a well defined a. The re¬ 
sult of such measurement in a Type-B device is shown 
in Fig. 3(b) for d = 400 nm between Fi and F 2 . Al¬ 
though the measured NLSV signal [Fig. 3(a)] is smaller 
than in Type-A devices, possibly due to a better Al/YIG 
interface, R^i exhibits a cos(2a) behavior with a maxi¬ 
mum (minimum) for a = 0 (a = 7r/2), consistent with 
Eq. (1). However, the maximum change (modulation) of 
the signal SRs=Rni{a = 0) — i?ni(cr = tt/2)) is only 12% 
of the total spin signal Rsv, which is at odds with the 
large spin-mixing conductance estimated from Fig. 2(b). 
From anistropic magnetoresistance measurements we ex¬ 
clude the possibility of any rotation of the magnetization 
of the injector and detector as the cause for the observed 
modulation in the NLSV signal (see Supplemental Mate¬ 
rial [19], Sec. III-B). 

Using the 3D-FEM we calculated the angular depen¬ 
dence of Rsv for various values of Gr where the percent¬ 
age modulation SRs/Rsv is plotted as a function of a, 
as shown in Fig. 3(c). The Gr value of 1 x 10^^ 
extracted from the NLSV signal modulation experiment 
is one order of magnitude less than reported elsewhere 
[4] . This can be possibly caused by the presence of disor¬ 
dered Al/YIG interface with r.m.s. roughness of 0.8 nm 
(as measured by AFM), which is close to the magnetic co¬ 


herence volume \fVc — 1.3 nm [6] of the YIG. This length 
scale determines the effective width of the Al/YIG inter¬ 
face and also the extent to which spin current from the 
A1 is felt by the YIG magnetization [6, 24]. Furthermore, 
the fact that there exists a finite spin-mixing when a = 0, 
as discussed above, can also explain the observed small 
modulation. It is important to note that in our experi¬ 
ments the non-equilibrium spin accumulation induced by 
electrical spin injection into A1 has a spin-polarization 
strictly along the direction of the magnetization of Fi, 
which lies along the y axis. In the measurement results 
shown in Figs. 1(b) and 2(b) the magnetization of the 
F 2 is always kept either parallel or antiparallel to the de¬ 
tector Fi. This ensures that it is only the y component 
of the spin accumulation that is measured in our exper¬ 
iments as it is insensitive to other two spin-polarization 


x-component z-component 



FIG. 4. (Color online) Calculated NLSV signals showing the 
(a) x-component and (b) 2 -component of the NLSV signal 
Rni in the parallel (red) and antiparallel (blue) magnetiza¬ 
tion configurations of the injector and detector ferromagnetic 
contacts for Gr = 1 x 10^® and Gi = O.IG,-. Even if 

the injected spin accumulation is polarized along the magne¬ 
tization direction of the injecting electrode Ei, its interaction 
with the magnons via the spin-mixing conductance induces 
these spin accumulation components. 
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directions. It is however possible that the interaction 
of the initially injected spin accumulation with the YIG 
magnetization, via to induce a finite NLSV signal 
with components polarized along the x- and z-directions. 

Figure 4 shows the angular dependence of the x— and 
z— component of the NLSV signal as calculated using our 
3D-FEM. While the z component exhibits a sin(a) de¬ 
pendence, the X component shows a sin(2a) dependence 
which is consistent with Eq. (1). The size of the mod¬ 
ulation is determined by G^ for the x— component and 
by Gi for the z— component. In a collinear measure¬ 
ment configuration these transverse spin accumulation 
components can induce local magnetization dynamics by 
exerting a spin transfer torque to the YIG. Separately 
measuring these spin accumulation using ferromagnetic 
contacts magnetized along the x and z directions can be 
an alternative way to extract G^^. 

In summary, we studied spin injection and relaxation 
at the Al/YIG interface in Ni 8 oEe 2 o/Al lateral spin valves 
fabricated on YIG. The samples on the YIG substrate 
yield NLSV signals that are two to three times lower than 
those grown on standard Si 02 substrates, indicating spin- 
current absorption by the magnetic YIG substrate. We 
also observed a small but clear modulation of the mea¬ 
sured NLSV signal as a function of the angle between 
the spin accumulation and magnetization of the YIG. 
The presence of a disordered Al/YIG interface combined 
with a spin-flip (sink) process due to thermal magnons or 
interface spin-orbit effects can be accounted for this small 
modulation. Using finite element magnetoelectronic cir¬ 
cuit theory as well as additional control experiments, we 
establish the concept of collinear (effective) spin mixing 
conductance due to the thermal magnons in the YIG. 
Our result therefore calls for the inclusion of this term in 
the analysis of spintronic and spin caloritronic phenom¬ 
ena observed in metal/YIG bilayer systems. 
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technical assistance. This work is part of the research 
program of the Eoundation for Fundamental Research on 
Matter (FOM) and is supported by NanoLab NL, EU- 
FET Grant InSpin 612759 and the Zernike Institute for 
Advanced Materials. 
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SUPPLEMENTAL MATERIAL 

I. Derivation for the effective spin 
relaxation length in the collinear case 


The spin accumulation fig, with polarization parallel to 
the magnetization direction of Fi (see Fig. S5), injected 
in the A1 is governed by the Valet-Fert spin diffusion 
equation [25] [9^ + dy + 9^] fig = which can be 

re-arranged to give 

dlfig = Hg/X% - d^^Hg. (5) 


Here we assume that, for a homogeneous system, the 
spin current along the y-direction is zero. As discussed 
in the main text, when the YIG magnetization direction 
TO |j fig the spin current at the Al/YIG interface, in 
the z-direction, is governed by the spin sink term Gg in 
Eq. 1 of the main text. Applying spin current continuity 
condition at the Al/YIG interface we find that 

^d^fig = Ggfig ( 6 ) 


where (Tat is the conductivity of the normal metal. Now 
after re-arranging Eq. (6) to obtain dzfig, differentiat¬ 
ing it once and using dzfig = fig/tAi, where Iai is the 
thickness of the Al, we obtain 


dzfJ^s 


‘IG gfig 

CTNtAl 


(7) 


Substituting Eq. (7) into Eq. (5) we obtain a modi¬ 
fied VF-spin diffusion equation that contains two length 
scales 


— Ms /-\- ‘IGj-fLg/OAitAl-, 


_ Ms I Ms 


(8a) 

(8b) 


where we defined a new length scale = 2GgjaNtAi 
that, together with the Aat, re-defines an effective spin 
relaxation length A“g = X~^^ + X~^. This effective spin 
relaxation length in the Al channel is weighted by the 
spin-mixing conductance Gg of the Al/YIG interface. 
The modulation of the NLSV signal observed in our mea¬ 
surements is hence determined by the interplay between 
these two length scales, Xn and A^. While the first quan¬ 
tifies the effective spin-conductance of the Al channel 
(Gat = (JpfApf/X ai) over the spin relaxation length, the 
second is a measure of the quality of the Al/YIG inter¬ 
face and is set by Gg. For the devices investigated in 
this work, using = tAiWAi with the width of the Al 
channel wai = lOOnm and gai = 2 x lO^S/m, we ob¬ 
tain A'^Gn ~ 6 X which is close to the Gg 

obtained in our experiments. This highlights the impor¬ 
tance of spin-relaxation induced by the thermal motion 
of the YIG magnetization, as discussed in the main text. 


Geometrical enhancement of the modulation can be ob¬ 
tained by reducing Iai, as shown in Fig. S5(d), thereby 
maximizing spin-absorption at the Al/YIG interface [16]. 


II. Three dimensional (3D) spin 
transport model 

Here we describe the our 3D spin transport model used 
to analyze our data. It is similar to that described in 
Ref. 20 for collinear spin transport with the possibility of 
studying spin-relaxation effects (i) due to the spin-mixing 
conductance at the Al/YIG interface as well as (ii) Hanle 
spin-precession due to the in-plane magnetic field [see 
Sec. HI below for detail]. The charge current /“(r) and 
spin current for j^{r} (where a G x,y,z), are re¬ 

lated to the charge and spin potentials fig as 


CT aFG\ fVflA 
\oiFG a ) VV/xJ 


(9) 


where a is the bulk conductivity and is the bulk spin 
polarization of the conductivity. The device geometry we 
model is shown in Fig. S5(a), showing schematic source- 
drain configurations as well as voltage contacts. We im¬ 
pose charge flux at contact 1 and drain it at 2. The 
nonlocal voltage, due to spin diffusion, is obtained by 
taking the difference between the surface integrated fi,. 
at contacts 3 and 4 , both for the parallel (P) or antipar¬ 
allel (AP) magnetization configurations. To solve Eq. (9), 
we use conservation laws for charge (V • (r) = 0) and 

spin current (V • jg = (1 — ap)Gfis/X^) with spin pre¬ 
cession due to the in-plane applied field also included in 
the model. By defining an angle a between fig and the 
YIG magnetization to and allowing for a boundary spin 
current at the Al/YIG interface using Eq. 1 of the main 
text, we can study the transport of spins in NLSV de¬ 
vices and their interaction with the YIG magnetization. 
The material parameters for the model, a, aF and A^ are 
taken from Ref. [23]. Our modeling procedure involves, 
first, fitting of the measured NLSV signal on a Si02 sub¬ 
strate by varying uf and using Aat = 350nm. Next, we 
aim to find Gg of the Al/YIG interface that properly 
quantifies spin transport properties of the YIG sample. 
Figure S5(b) shows the dependence of the NLSV signal 
on Gg. As expected, when Gg very low, the NLSV signal 
is not affected by the presence of the YIG as spins are 
not lost to the substrate. For Gg ~ 5 x we 

obtain the experimentally measured NLSV signal (shown 
in red dashed line). For even larger Gg values, the effect 
is maximum with the NLSV signal falling by almost one 
order of magnitude. It is important to remember that 
the value of Gg that is extracted here is a simple mea¬ 
sure of spin-flip processes at the Al/YIG interface due to 
thermal fluctuation of the YIG magnetization or disorder 
induced effects. At the temperatures of our experiment it 
is difficult to distinguish which one of the two processes 
is dominant. 
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FIG. 5. (a) Geometry of the modeled device showing the measurement configuration with a 3D profile and the y-component 
of the spin accumulation, (b) The dependence of the NLSV signal on the effective (collinear) spin mixing conductance Ga- To 
reproduce the experimentally observed decrease in the spin signal from Si02 to the YIG substrate, an effective spin mixing 
conductance of Gs = 5 x 10^® is required, (c) The dependence of the NLSV signal on the angle between m and fia 

for Gs = 5 X (d) The dependence of the spin signal modulation amplitude on the thickness of the A1 channel 

signifying the interplay between the spin-mixing conductance and the spin-conductance in the A1 channel. 


For the angular dependent simulation we only vary the 
angle a between fig and rh while keeping all other param¬ 
eters constant (such as ap, Ajv = 5 x 10^^ 
and Gr = 8x 10^^ As shown in Fig. S5(b) our 

simulation as described above reproduces the cos^(a) de¬ 
pendence observed in our experiments as well as by Vil- 
lamor et al. [16]. 

For the extracted values of Gr from our analysis, the 
experimentally observed modulation of the NLSV signal 
by the rotating magnetization direction of the YIG is 
small. Possible ways to enhance the modulation are to 
1 ) maximize the spin-mixing conductance via controlled 
interface engineering of the Al/YIG interface or 2) reduce 
the thickness of the spin transport channel. In the latter, 
for a fixed Gr, the effect of decreasing the thickness of 
the spin transport channel is to effectively reduce the spin 
conductance Gn along the channel thereby maximizing 
the spin current through the Al/YIG interface. Fig¬ 
ure S5(c) shows the thickness dependence of the modula¬ 
tion of the spin signal SRs = Rs^a = 0^) — Rsia = 90°) 


normalized by Rg as a function of the thickness Iai , with 
the inset showing that for the P and AP configurations. 
As the thickness of the Al channel increases the spin cur¬ 
rent absorption at the Al/YIG interface decreases or vice 
versa. 


III. Investigation of possible alternative 
explanations for the observed modnlation 

It can be argued that the experimentally observed 
modulation of the NLSV signal can be fully explained 
by (i) the Hanle spin-precession and/or (ii) the rotation 
of the magnetizations of the injector/detector electrodes 
due to the 5 mT in-plane magnetic field. Below, we show 
that even the combined effect of both mechanisms is too 
small to explain the experimentally observed modulation 
of the NLSV signal. 
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FIG. 6. (a) Modulation of the NLSV when only considering the Hanle effect due to the in-plane magnetic field in the P 
(dashed lines) and AP (solid lines) at 5 mT (red), 50 mT (blue) and 100 mT (black), see text for more details, (b) Anisotropic 
magnetoresistance (AMR) measurement for the injector (left) and detector (right) ferromagnets at two different magnetic fields. 
The insets show the full-scale plot of the measurements at 5mT. 


A. Hanle spin-precession induced modulation of the 
NLSV signal 

Spins precessing around an in-plane magnetic field 
B would acquire an average spin precession angle of 
(j) = loltdi where wl = g^sB/h is the Larmor precession 
frequency, rp = L? j 2D ^ = 25 ps is the average diffusion 
time an electron takes to traverse the distance L between 
the injector and the detector and D(. = 0.005m^/s is the 
diffusion coefficient [26]. For an applied field of 5 mT 
and L=500 nm, we obtain (f> = 1.25°, giving us a max¬ 
imum contribution of 1 — cost/) =~0.02% [see Eq. (10)] 
to the experimentally observed signal (compared to the 
~12% in Fig. 3(b) of the main text). This is expected be¬ 
cause the spin-precession frequency (~8 ns) at such 
magnetic fields is three orders of magnitude slower than 
td- 

This simple estimate is further supported by our 3D fi¬ 
nite element model as we show next. Figure S6(a) shows 
the angle dependence of the nonlocal signal due to an 
in-plane magnetic field when we only consider the Hanle 


effect both for the AP (solid lines) and P (dashed lines) 
configurations at three different magnetic field values of 
5 mT (red), 50 mT (blue) and 100 mT (black). The 
maximum modulation of the NLSV signal that the Hanle 
effect presents is only 0.001% at the measurement field 
of 5 mT and only become relevant at high fields. There¬ 
fore, the Hanle effect alone can not explain the results 
presented in the main text. 

B. Magnetization rotation induced modulation of 
the NLSV signal 

The in-plane rotation of the sample under an applied 
magnetic field of 5 mT might induce rotations in the ma- 
gentization of the injector/detector electrodes. In such a 
case, a relative angle 9r between the magnetization direc¬ 
tion of the injector and detector electrodes would result 
in a modulation of the NLSV signal given by 

5Rnl _ Rnl{dr = 0 ) — Rnl{9r) _ ^11 a I 

Rnlidr = 0) " = 0) “ ' -COSy.l, 

( 10 ) 
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with +(—) corresponding to the P (AP) configuration. 
Using Eq. (10), we find that a relative angle 9r — 28° 
between the magnetization directions of the injector and 
detector is required in order to explain the experimentally 
observed modulation. To determine the field induced in¬ 
plane rotation of the magnetization by the applied mag¬ 
netic field, we carried out angle dependent anisotropic 
magnetoresistance (AMR) measurements both for the in¬ 
jector and detector electrodes, using a new set of devices 
with identical dimensions. The AMR measurements were 
repeated for different magnetic field strengths, at 5 mT 
and at higher magnetic fields of 100 mT and 300 mT. 

Figure S6(b) and (c) show the two-probe AMR mea¬ 
surement of the injector and detector electrodes, respec¬ 
tively, at two different magnetic fields. For the injector 
electrode in Fig. S6(b), at an applied field of 100 mT 
(red line), an AMR response AR = Ry — R± = 0.6 U 
is observed, where i?||(i?x) is the resistance of the ferro- 
magnet when the angle between the applied field and the 
easy axis is 0 = 0° (0 = 90°). For the same electrode, at 
an applied field of 5 mT (blue line, see also the inset), 
the AMR response is only 0.025 fl. Now, by comparing 
these two measurements we conclude that the effect of 
the 5 mT field would be to rotate the magnetization of 
this electrode by a maximum angle 0i = 15° from the 
easy axis. A similar analysis for the detector electrode, 
using the AMR responses of 2 U (at 300 mT) and 0.025 
U (at 5 mT) in Fig. S6(c), yields a maximum magneti¬ 
zation rotation 02 = 10°. Relevant here is the net rel¬ 
ative magnetization rotation between the two electrodes 
0r = 01 — 02 = 5° and, using Eq. (10), we conclude that 
it would only cause a modulation of 0.4 %, which is much 
smaller than the 12% observed in our experiments. Our 
analysis based on the AMR effect is equivalent to that in 
Ref. 16 where magneto-optical Kerr effect measurements 
were used to exclude a possible in-plane magnetization 
rotation as the origin for the observed modulation in the 
nonlocal spin valve signal [16]. 

To summarize this section, the Hanle effect and the 
magnetization rotation induced by the in-plane magnetic 
field neither separately nor when combined are sufficient 
to explain the experimentally observed modulation. Only 
after including the effect of the spin-mixing interaction 
via that it is possible to reproduce the modulation 
observed in the experiments. 


* e-mail:f.k. dejene@gmail.com 
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